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We study the diffraction of quantum degenerate fermionic atoms off of quantized light fields in 
an optical cavity. We compare the case of a linear cavity with standing wave modes to that of a 
ring cavity with two counter-propagating traveling wave modes. It is found that the dynamics of 
the atoms strongly depends on the quantization procedure for the cavity field. For standing waves, 
no correlations develop between the cavity field and the atoms. Consequently, standing wave Fock 
states yield the same results as a classical standing wave field while coherent states give rise to a 
collapse and revivals in the scattering of the atoms. In contrast, for traveling waves the scattering 
results in quantum entanglement of the radiation field and the atoms. This leads to a collapse and 
revival of the scattering probability even for Fock states. The Pauli Exclusion Principle manifests 
itself as an additional dephasing of the scattering probability. 
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I. INTRODUCTION 

The past few decades have witnessed considerable 
progress in the cooling of atomic vapors to extremely low 
temperatures, culminating in the achievement of Bose- 
Einstein condensation in dilute alkali gases LJ- LjI ■ More 
recently, quantum degenerate Fermi gases with temper- 
atures as low as 0.017>, where Tp is the Fermi temper- 
ature, have been achieved by several groups 0, 01 El 
Throughout these developments the interaction of light 
with atoms has been central to the cooling, trapping, and 
imaging of atoms, as well as in the coherent manipula- 
tion of their center-of-mass motion. For example, the 
Bragg scattering of atomic matter waves by off-resonant 
optical fields can be used to create linear atom optical 
elements for use in atom interferometers Q, and the 
interaction of atomic condensates with light has led to 
the realization of matter-wave superradiance @ and of 
matter- wave parametric amplifiers 0, ITol Hlj . In another 
application, the ability of optical fields to create custom 
trapping potentials has permitted the study of condensed 
matter problems such as e.g. the Mott-Insulator transi- 
tion [lj) EU EH • Although all experiments to date have 
involved classical optical fields, there is considerable in- 
terest in carrying out future work in high-Q optical cav- 
ities, where the quantum nature of the electromagnetic 
field becomes important. Theoretical work along these 
lines has so far been restricted to the case of bosonic 
atoms, see e.g. Ref. [l5j . while the diffraction of fermions 
by an optical field was discussed in Ref. |16|. but in 
an analysis restricted to the case of classical fields. In 
this paper we extend this work to discuss the diffraction 
of quantum-degenerate fermionic matter-wave fields by 
quantized light fields. 

We consider a zero-temperature beam of fermionic two- 
level atoms traversing an optical cavity supporting an 
off-resonant standing wave light field of momentum q. 
The atoms undergo virtual transitions to their excited 
electronic state, resulting in a center-of-mass momentum 



recoil of 2q . Alternatively, one can view this process as 
diffraction of the atoms off of the intensity grating formed 
by the cavity field. 

The normal modes in terms of which the electromag- 
netic field is quantized are determined by the boundary 
conditions of the cavity. In a linear cavity with per- 
fectly reflecting mirrors we have standing wave mode 
functions. In a ring cavity, the light field has to fulfill 
periodic boundary conditions and this results in running- 
wave mode functions. Two counter-propagating travel- 
ing wave modes of equal frequency can be superposed to 
yield a stationary standing wave field. 

Under most circumstances it is a question of mathe- 
matical convenience which mode functions are used for 
the description of the field. Physically, however, the two 
cases are not the same and for fields containing only a few 
photons, the two quantization procedures yield different 
results. In particular, the difference in atomic scattering 
produced in these two situations has been discussed for 
single atoms diffracted by a coherent light field ^t|- ^ 
was shown to depend critically on the quantization pro- 
cedure. This difference can be understood in terms of 
which-way information for the scattering process. For 
standing wave modes, the state of the light field con- 
tains no information about the momentum transfer to 
the atom. More specifically, the number of photons is a 
constant of motion and as a result the equations of mo- 
tion for the atomic center of mass decouple from that of 
the light field. In the case of two counter-propagating 
traveling wave modes however, the number of photons in 
each mode does change and the change in the number 
of quanta is a direct measure of the momentum transfer 
to the atoms. In this paper we extend those results to 
compare the diffraction of a quantum-degenerate Fermi 
gas by these fields, both in the Raman-Nath and Bragg 
regimes. 

In the Raman-Nath regime, which is characteristic of 
situations where the kinetic energy of the atoms can be 
neglected, the individual atomic dynamics for a stand- 
ing wave light field are formally identical to the case of 
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a classical light field ^^|. The atoms scatter into suc- 
cessive diffraction orders separated by twice the photon 
momentum q, up to the point where energy-momentum 
conservation becomes important and the Raman-Nath 
approximation ceases to hold. The formal equivalence of 
the scattering off of a standing wave field to the scat- 
tering off of a classical light field is due to the fact that 
the equations of motion for the atoms effectively decou- 
ple both from each other and from the light field. This 
must be contrasted to the case of running waves, where 
the number operators for the two modes are not con- 
stants of motion. This leads to an infinite hierarchy of 
coupled equations for the atomic and optical field oper- 
ators, with higher-order correlation functions playing a 
crucial role in the dynamics of first-order atomic corre- 
lation functions. It is then necessary to introduce some 
approximate truncation scheme, a procedure that we dis- 
cuss in detail and compare with exact numerical results 
for small atom numbers. 

In the Bragg regime, energy-momentum conservation 
reduces the single-atom diffraction problem to a two- 
mode situation, the atoms undergoing Bragg oscillations 
between their initial momentum states, pi, and final mo- 
mentum state pf — pi + 2q. The character of these os- 
cillations is the result of three separate and independent 
effects which correspond to whether one uses standing 
wave or traveling wave modes, whether the cavity field is 
in a Fock state or in a coherent state, and the momentum 
spread of the incident atomic beam. 

This paper is organized as follows: After formulating 
the specific model used in our analysis in section [n] wc 
discuss the case of traveling-wave light quantization in 
section IIIII We develop approximate equations for first 
and second-order correlation function appropriate for the 
Raman-Nath regime, and a Bloch vector picture useful to 
discuss Bragg diffraction. Specifically, that picture yields 
a semiclassical model that provides some intuitive under- 
standing of the atomic dynamics. The case of standing- 
wave quantization is discussed in Section llVl Section W\ 
gives a summary and conclusion. 



II. MODEL 

We consider an ultracold beam of identical two-level 
fcrmionic atoms propagating across a high-Q optical cav- 
ity, see Fig. ^ Their initial momentum distribution is a 
Fermi sea at T = 0, but shifted in momentum space by 
the mean momentum p and with Fermi momentum fcp 
assumed to be much less than q, the photon momentum. 
This is a realistic approximation, since for a degenerate 
Fermi gas of density n « 10 17 m -3 the Fermi momen- 
tum is kp ps 10 6 m _1 while for a photon of wavelength 
A = 500nm one has a momentum of q w 10 7 m _1 . In the 
following we neglect atomic collisions - a good approxi- 
mation at low temperatures, since the s-wave scattering 
length is zero for identical fermions - as well as cavity 
losses. We also assume that the optical frequency lu is 
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FIG. 1: (a) Schematic of a scattering of an atom of initial 
momentum p via two photon transitions with photons of mo- 
menta q and —q. (b) Initial momentum distribution N(p) of 
the atoms for scattering in the Raman-Nath regime and (c) 
for the Bragg regime. 



sufficiently detuned from the atomic transition frequency 
wo that the upper electronic level can be adiabatically 
eliminated. Finally, we consider a situation where that 
atomic momentum mv± transverse to the cavity field is 
large enough so that it can be treated classically. Time, 
t, can then be parameterized in terms of the transverse 
distance x by t = x/v±. 



III. 



RUNNING WAVES 



For running-wave quantization, the Hamiltonian de- 
scribing our system is {h = 1) 
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uj(a\a q 
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where, c k and c k are the annihilation and creation op- 
erators for a fermionic atom of momentum k, a q and a q 
arc the annihilation and creation operators for a photon 
of momentum q, E k = k 2 /2M is the kinetic energy of an 
atom of momentum fc, g = flj^/A is the coupling energy 
of the atoms and the light field, £Ir is the vacuum Rabi 
frequency, and A = uj — ujq is the atom-light detuning. 
The initial state of the atoms-field system is 



h/>(o)> rw = i^>i*_ fl ) n 4io). 



(2) 



|fc|<fci 



where the field states \4>± q ) are taken to be either Fock 
states \N± q ) or coherent states \a± q ). 

A. Raman-Nath regime 

The Raman-Nath regime of atomic diffraction is char- 
acteristic of situations where the kinetic energy of the 
atoms plays a negligible role in comparison with the in- 
teraction energy, i.e. E 2q <C g\J 'N q N^ q , where the recoil 
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energy Ei q = 2q 2 /M is a measure for the typical kinetic 
energies involved. In practice, this amounts to assuming 
that the atoms have an infinite mass, and as such, ne- 
glects the effects of the quadratic dispersion relation of 
the atoms. 

The most straightforward way to solve this problem 
proceeds by integrating the Schrodinger equation corre- 
sponding to Hamiltonian Q for the initial conditions J2J, 
From which we can obtain the probability 

p p (t) = m)\4>cp\m)- (3) 

for an atom being scattered to a state of momentum p. 
However, the dimension of the Hilbert space grows ex- 
ponentially as Dim Raman „ Nath = (2n d + l) N »(N p + 1) 
where rid is the number of diffraction orders considered, 
iV a is the number of atoms, and N p the total number of 
photons. Hence, a direct integration of the Schrodinger 
equation is only possible for rather small atom and pho- 
ton numbers. 

Figure[2Ia) shows the result of an exact solution of the 
Schrodinger equation for N a = 2 and the light held in 
a Fock state with three photons per mode initially. In 
this example, the recoil energy is Ei q = g and the ini- 
tial momentum of the atoms is pi — ±0.1q. Such a high 
recoil energy was chosen to limit the number of diffrac- 
tion orders that are significantly populated before energy- 
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where the 6's are Kronecker-deltas. 

The simplest factorization scheme consists in merely 
factorizing second-order correlation functions of the type 
(a^ciqzC^Cfo) that appear on the right-hand side of these 
equations into products of first-order correlation func- 
tions, for instance, (a^a^cl^) « (<4i a <s)(4i c *2) ■ 
In doing so we neglect correlations that may build up 
between the atoms and the light field as well as higher- 
order correlations of both the atoms and the light field. 
This corresponds to a truncation of the BBGKY-type hi- 
erarchy of the equations of motion for the higher-order 
moments of the particle-hole operators after the first or- 
der, see e.g. This reduces the infinite hierarchy of 
equations 1 llol) to a closed set of c-number equations that 
grows only quadratically with the number of momentum 
states that have to be taken into account. 

The result of the numerical integration of these equa- 
tions of motion, using the same parameters as previously 



momentum conservation inhibits further diffraction, i.e. 
before exiting the Raman-Nath regime. 

The resulting dynamics resembles qualitatively the 
single-atom case, see e.g. [l^. For short times the prob- 
ability for finding an atom in the mth order mode is well 
described by ~ J 2 n (2gt) where J m is the m-th Bessel 
function. For longer interaction times, higher scattering 
orders are suppressed due to energy-momentum conser- 
vation, as expected. We note that since the difference in 
kinetic energies of the two atoms is small compared to 
all other relevant energies, we do not observe any effect 
of "inhomogeneous broadening." 

For comparison, the results for initial coherent states 
with mean photon numbers N q — N_ q — 3 are shown 
in Fig. |2fb), the atomic parameters being the same 
as before. We now observe a decay of the oscilla- 
tions of the scattering probabilities after a time t ~ 
(2tt I g)(N q N -q)^ 1 / 2 , which corresponds to a complete 
dephasing of the contributions of the different photon 
numbers to the diffraction pattern. 

In order to proceed past the few-atom problem, we now 
concentrate on single-particle properties, introducing a 
BBGKY-type truncation scheme to factorize higher- 
order correlation functions of the matter-wave field. 
From the Hamiltonian {T|l. the equations of motion for 
the atomic first order correlations (c£ Cfe 2 ) and (a qi a q2 ), 
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for ease of comparison, is shown in Fig[2c,d) for the 
cases of a Fock state and a coherent state of the field, 
respectively. 

An obvious weakness of the simple truncation scheme 
is that it predicts the absence of scattering for the case 
of Fock states, in stark contrast to the exact solution. 
This follows from the absence of initial coherence in ei- 
ther the light field or the atoms, leading to the scattering 
term in <0J being identically zero. Stated differently, the 
reason for the absence of diffraction is that the phase of 
a Fock state is completely undetermined, hence there is 
no established relative phase between the two counter- 
propagating fields, and no light intensity grating. Since 
in this factorization scheme the atom is effectively as- 
sumed to probe only first-order moments of the light field, 
that is, its intensity pattern, diffraction is absent at this 
level of approximation. 

The situation is different for a coherent state light field. 
In this case, there is a well-established phase relation- 



(E k2 - E kl )(c{ iCk2 ) 



g (al q a q Uck^q - c{ i+2q c kl 



+ 9 (4^2+29 ~ 4 1 -2g C *2) ) , (4) 



9j2( 5 q2, q a f qi a- q c{_ q c k+q + 5 q2 ^ q a\ iaq cl +q c k - q + 5 qi ^ q a\a q2 cl_ q c k+q + 8 qi , q al q a q2 c{ +q c k - q ^ , (5) 




FIG. 2: Scattering probability P p (t) for two atoms scattering off of a running- wave light field in the Raman-Nath regime. 
Figures (a,c,e) are for a Fock state of the light field and figures (b,d,f) for a coherent state. Figures (a,b) show the exact 
solution of the Schrodinger equation, figures (c,d) for the first-order equations and figures (e,f) show the results for the second- 
order equations. In all the calculations the recoil energy is E-2 q = g and the Fermi momentum is &f = 0.1. Time is in units of 
g^ 1 and momentum in units of the photon momentum q. 



5 



ship between the two modes. This results in an inten- 
sity grating from which the atoms can be diffracted. As 
time goes on, this results in the generation of atomic co- 
herence, (cj^Cfcj,) ^ 0, and the resulting density grating 
formed by the atoms acts back on the light field. In 
some loose sense, the lowest order factorization scheme 
consists in treating the system classically since it neglects 
all quantum fluctuations in the atomic and optical fields. 
It is not surprising that this approach should fail for a 
very non-classical field state such as a Fock state, and 
be much better for a quasi-classical field. Note however 
that while for short enough times the scattering closely 
resembles the exact results, this is no longer the case for 
long times, a consequence of the build-up of quantum 
correlations between the optical and matter-wave fields. 
Even after one oscillation differences arise. 

We note that for our specific initial conditions, the 
fully factorized equations for the light field can be triv- 
ially integrated, showing that the first-order moments of 
the light field are constants of motion. Inserting these 
constants in the atomic equations of motion shows that 



at this level, the scattering becomes formally equivalent 
to the scattering of atoms by a classical standing wave 
light field with intensity g(a q a- q ). This is further dis- 
cussed in the following section. 

The equations of motion (1 1I5H . suggest that an im- 
proved factorization scheme would retain the lowest or- 
der correlations between light field and atoms. In or- 
der to do so, we supplement the equations of motion for 
(c' k Ck 2 ) and (a qi a q2 ) by equations of motion for the cross- 
correlations (aj 1 a 9 2 c i 1 c fc 2 ) an d the second order correla- 
tions of the lightficld and the atoms. This should remedy 
the major flaw of the first-order calculation, namely its 
inability to predict atomic scattering for a light field in a 
Fock state. 

The equations for the lowest-order atom-field correla- 
tion functions involve third-order correlations of the form 
(4i tt » fl l 94 c L c fc2> and (4i a <&4i c *2 c L c *4)- We trun- 
cate the resulting hierarchy of equations of motion by 
introducing the factorization scheme 
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and similarly for (q^ a Q „cl c k 2 cj„ c k±) with a q s replaced 
by Cfc's and vice versa |20j. The last term of this equa- 
tion accounts for the case where all first-order correlation 
functions are uncorrelated. 

We estimate the accuracy of the factorization 
scheme by calculating {a\ 1 a q 2 a \* i a qi c \ 1 c k 2 ) and 
(a qi a q2 c' ki Cfa c£ Cf-i ) as well as their respective fac- 
torized values, Eq. (jSJ) using the exact solution of the 
Schrodinger equation. As an example Fig. |3| shows the 
results for {a}_ q a q a q a- q c kp Ck F ) and Fig. 0] shows the 

results for (a q a q c kF Ck F c k Ck P ) for the parameters of Fig. 
121 This shows that the factorization scheme reproduces 
at least qualitatively the main features of the third-order 
correlation functions for both coherent states and Fock 
states. 

Despite its apparent success, we must keep in mind 
that this factorization scheme suffers from two major 
flaws. First, the small deviations of the factorized val- 
ues from the exact values will accumulate in the course 
of time, leading to increasing discrepancies between the 
approximate and exact results. More critical perhaps, 
this scheme violates important relations that the exact 
operators have to obey. For example, CL q a q c s kF Ck F c s kF Ck F 
is a positive self-adjoint operator, with positive and real 
expectation values, but the factorized approximation can 
take on negative values. These flaws eventually result in 
non-physical behavior such as illustrated in Fig. H[e,f) 
where the probabilities (c k Ck) take on negative values. 



We have not found a factorization scheme that avoids 
non-physical behavior of that kind at all times, and con- 
jecture that the factorization of higher-order moments in 
lower-order moments necessarily leads to such inconsis- 
tencies. 

The results of the factorization scheme @ are shown 
in Fig. 0e) for a Fock state and in Fig. (2ff) for a co- 
herent state of the light field. While a Fock state now 
leads to atomic diffraction, as should be the case, it is 
characterized by non-physical negative probabilities al- 
ready for short times. This is clear evidence that higher 
order correlations play an essential role. This is in con- 
trast with the situation for a coherent state, where we 
achieve good agreement with the exact results for times 
up to ~ (27r / 'g)(N ~ q N ^q) -1 / 2 , indicating that the first 
and second-order correlations are the most important. 

A quantitative measure of the degree of entanglement 
between the atoms and the light field is given by the 
second-order cross-correlation 

= (( a l a -1 C k~q C k+q) - ( a la-q)( C i-q C k+q)) , 

k 

which is equal to zero in the absence of entanglement. 
Figure El shows x(t) for both a Fock state and for a co- 
herent state light field, for the parameters of Figj2J Be- 
cause the light field and atoms are initially uncorrelated 
we have x(t = 0) = but cross-correlations then build 
up to become of the order of (N g N ^q) 1 / 2 N a . The fig- 
ure also shows the result of the factorization ansatz ||§J , 
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FIG. 3: Expectation value of a'_ q aqa,qa-qC' kp Ck F as calculated 
from the numerical solution of the Schrodinger equation for 
a Fock state and for a coherent state of the light field with 
the same parameters that were used in Fig. [5] o and o show 
the exact unfactorized value for Fock state and coherent state 

respectively; the broken line ( ) and the solid line ( — ) 

show the corresponding values for the Fock state and coherent 
state as obtained with the factorization scheme (JHJ- 




FIG. 4: Expectation value of a q a q c' kF Ck F c' k Ck F as calculated 
from the numerical solution of the Schrodinger equation for 
a Fock state and for a coherent state of the light field with 
the same parameters that were used above, o and o show 
the exact unfactorized value for Fock state and coherent state 

respectively. The broken line ( ) and the solid line ( — ) 

show the values obtained from the factorization scheme © 
for the two cases. 



showing the good agreement with the exact result for 
short enough times. 

B. Bragg regime 

In the Bragg regime, energy-momentum conservation 
restricts the scattering of the atoms to two diffraction 
orders, an initial mode of transverse momentum w — q 
and a final mode of momentum pf = pi + 2q « q. Classi- 
cally, the atoms are known to undergo Pendellosung os- 
cillations between these two modes. As such, the atoms 
can be thought of as two-state systems that are conve- 
niently described in terms of pseudo-spin operators 

&k — 2 ( C k+q C k+q - C \-q C k~q) 

st = is^ = i +q ^- q - (8) 

Introducing further the Schwinger representation of the 
light field by means of 

J" = 2 ( a l a 1 ~ a -q a -l) ' 

J+ = (J-)t=ota_„ (9) 
the Hamiltonian of the atoms-field system simplifies to 
H= J2 {6ukS* k +gJ + S k -+gJ-S+). (10) 

k£ [—k F ,k F ] 

In this representation, the eigenvalues m of J z corre- 
spond to the photon number difference m = (l/2)(iV ? — 



N- q ) between the two counter-propagating modes and 
the eigenvalues + of J 2 = J z2 + l/2(J+J- + J-J+) 
correspond to the total number of photons, j = l/2(N q + 
N- q ). Finally, Scut = E k+q - E k - q = 2kq/M is the 
frequency mismatch between the two momentum states 
accessible to the atom with initial momentum k — q. 

When compared to the Raman-Nath case, the dimen- 
sion of the Hilbert-space is now reduced to Dimeragg = 
(N p + 1)2^°, which allows us to consider larger atomic 
and photon numbers. The initial state of the atoms is 
now a T = Fermi sea shifted by —q in momentum. 
In the following, we evaluate the total number of atoms 
diffracted to states of momentum near +q, 

fc£ [ — k F , k F ] 

for a light field initially in a Fock state and in a coherent 
state. 



1. Fock state 

Figure EJa) shows P Pf (t) for a Fock state with N q — 
N- q = 6 photons and five atoms, with a recoil energy 
E2 q = 50g and the Fermi momentum is kp — O.lq. Fig- 
ure[7[a) show N sc (t) for the same parameters. 

For short times the atoms undergo Pendellosung oscil- 
lations between initial and final momentum states, with 
an amplitude that decreases as the atomic momentum 
is further detuned from the Bragg resonance condition. 
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FIG. 5: Cross-correlation \(t) between light field and atoms in 
the Raman-Nath regime as obtained from the exact solution 
of the Schrodinger equation for a Fock state (o) and a coherent 
state of the light field (o). Also shown are the values obtained 

from the factorization scheme © ( ) for a Fock-state light 

field and ( — ) for a coherent state. Same parameters as in 
Fig. 2. 



For longer times the oscillations of the individual atoms 
dephase, as expected from their different kinetic ener- 
gies. However, the dephasing is not as strong as would 
be the case for a system of independent particles, such 
as in the standing-wave case discussed later on, a clear 
manifestation of the collective nature of the system. 

In addition, the individual atomic oscillations undergo 
a decay that is intrinsically linked to the quantum cor- 
relations that build up between the light field and the 
atoms and is present even if we neglect dephasing (i.e. 
fejr = 0), as shown in FigEJb). For zero dephasing, the 
amplitude of the oscillations eventually revives to its ini- 
tial value. (The fact that the collapse of the oscillations 
and their subsequent revival resemble a beat phenomenon 
in the figure is an artifact from the comparatively small 
number of atoms and photons.) Combined with the in- 
homogeneous dephasing due to the width of the Fermi 
sea, this decay results in the total oscillation amplitude 
N sc {t) shown in FigEfa). 

We can gain a qualitative understanding of the col- 
lapse and revival from an analysis of the matrix elements 
of the operators J + and J~ , which give an estimate of 
the transition frequencies for the atoms from pi to pf. 
although for an initial Fock state the system starts in 
a state of definite to, it evolves over time into a linear 
superposition of to states. The matrix elements of J + 
and J~ between different m-states yield different Rabi 
frequencies and hence Bragg oscillation periods. Eigen- 
states of J z with eigenvalues to and m + 1 are coupled 
by the matrix element 



We can therefore estimate the collapse time, Td eC ay, by 
calculating the difference between the fastest and the 
slowest of these frequencies, the collapse time being 
roughly the time after which this frequency difference has 
produced a phase difference of 2tt. Under the assumption 
that all TO-states contribute equally to the dynamics we 
find 



T, 



decay 



27T.9- 1 



(i,l|J+|i,0)-(j,j|J+|i,i-l) 



VUTTJj-VZj' 



(13) 



This estimate gives satisfactory agreement with the ac- 
tual decay time for small j (it is within ~ 10% of the nu- 
merical result for our parameters), but breaks down for 
large j. We attribute this to the fact that the assump- 
tion that all to states are initially equally populated is 
unphysical for large j. 

The revival time of the Pendellosung oscillations can be 
evaluated in a similar fashion: The revivals occur when 
the Rabi frequencies for neighboring TO-states differ in 
phase by 2ir. This gives 



cvival 



(j,l\J+\j,0)-{j,2\J+\j,l) 

Vti+W -y/{j + 2)(j-\y 



(14) 



(j,m+ l\J + \j, to) = + TO + l)(j - to). 



(12) 



which goes to infinity when j — > oo. Note that this esti- 
mate is not limited to small j values since it does not rest 
on the assumption of equal populations of all to states. 

The collapse and revival times can be evaluated more 
quantitatively from the spectrum of the Hamiltonian, as 
illustrated in Fig. [H] While the eigenfrequencies cover 
the whole spectrum rather densely, the initial state of 
the system is well described as a superposition of just 
a few groups of eigenstates, and hence only a few nar- 
row bands of frequencies, which turn out to be almost 
equally spaced, significantly determine the atomic dy- 
namics. If these frequency bands were exactly evenly 
spaced the Pendellosung oscillations would be perfectly 
periodic. The variations in spacing and the widths of the 
various frequency bands lead to the more complicated 
dynamics. 

The width of the frequency bands, which can be traced 
back to the usual dephasing of the atoms due to their 
spread in kinetic energies, gives the ordinary decay of 
the density oscillations, while the variation in separations 
between bands is a measure of the inverse revival time. 
Figure shows this separation, obtained numerically for 
several photon numbers with and without dephasing. For 
comparison we also give the inverse revival time as de- 
termined from the matrix elements of J + as well as by 
a direct inspection of N sc (t). While the agreement be- 
tween the revival times determined from the spectrum of 
the Hamiltonian and from N sc (t) is good for all photon 
numbers, the agreement with the estimate based on the 
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FIG. 6: Scattering probability P PJ (t) for five fermions scattering off of a running-wave light field (a) and a standing-wave light 
field (b) in the Bragg regime. In both cases the light field is in a Fock state with N q — N- q = 6 and N = 12 respectively, 
E2 q — 50g and &f = O.lg. Time in units of g~ and momentum in units of q. 
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FIG. 7: Mean number of scattered atoms, N ac = ^ P Pf for 
Fock states of the light field and for different recoil energies, (a) 
Running wave light field with N q = iV-, = 6 photons in each 
of the two modes with recoil energy i?2 9 = 50g. (b) Same light 
field as in (a) but without dephasing for the atoms, Ei q = 0. 
(c) Standing wave light field with N = 6 photons with recoil 
energy Ei q — 50<? for 200 atoms. Time in units of g 1 - 



FIG. 8: Projection of the initial conditions onto the spectrum 
of the Hamiltonian llOU for N q = N- q = 6 photons, five atoms, 
a recoil energy of E2 q = 50g and a Fermi momentum of k,F = 
O.lq. Eigenfrequencies in units of g. 



matrix elements of J + improves for large photon num- 
bers. 



2. Coherent state 

The case of a coherent state is readily obtained by 
averaging the Fock state results over a Poissonian photon 
distribution. The results for P Pj {t) and N sc {t) are given 
in Figs. HUf a') and II If a), respectively. The oscillations 



of the mean number of scattered atoms as well as the 
Pendellosung oscillations of the individual atoms decay in 
a time < g^ 1 . In addition to the effects discussed in the 
previous section, we now have an additional dephasing 
due to the photon statistics of the coherent states. These 
independent dephasing processes are normally associated 
with non-commensurate decay rates and revival times, 
hence there are no revivals in this case. 
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where we have introduced the total atomic spin operators 
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FIG. 9: Inverse revival time 27r/T rev ivai of the oscillations as 
determined from the spectrum of the Hamiltonian as a func- 
tion of the photon number for Ei q — 50g and &f = 0.1g(0)- 
For comparison we also show the results without dephasing, 
kF — (o) and the values obtained from comparing matrix 
elements of J+ , i.e. + Tjj - V '(j + 2)(j-l) (□). The 

inverse revival times as determined directly from simulations 
for N sc (t) like the one shown in Fig. [7Tb) for no dephasing 
are also given (V). Frequencies are in units of g. 



C. Bloch vector model 

If we wish to consider greater atom numbers we can 
introduce a factorization scheme such as we used in the 
Raman-Nath regime. Here, the fact that the individ- 
ual atoms are two-level systems in momentum space sug- 
gests that we use instead an approach analogous to the 
use of the Bloch vector in conventional quantum op- 
tics. We proceed by introducing the pseudo-spin vector 
Sk = S%y + S^z where x, y, and z are unit vectors 

along the x, y, and z directions in the abstract Bloch 
vector space and 



Equations i|17f) are exact within the two-state approxima- 
tion of Bragg scattering. We can obtain a semiclassical 
picture by factorizing expectation values of products of 
atomic and field operators, e.g. (S m J n ) = (S m )(J n ). 

For atoms with initial momenta centered around — q 
we have (5|(0)) = -1/2, (S*(0)) = (S*(0)) = for 
all k, so that the individual atomic Bloch vectors point 
to the south pole. Likewise, for a field in a Fock state 
we have that (J x (0)) = {J v (0)) = 0, so that J points 
along the z-axis, too. From the Bloch equations i|17|) 
it is then immediately apparent that there is no atomic 
scattering in the semiclassical description, consistently 
with the previous discussion. 

For a coherent state, on the other hand, J is not 
parallel to the z-axis. For our choice of phase and for 
N q = N_ q , it points instead along the x-direction. The 
phase relationship between the two counter-propagating 
coherent states leads to an intensity grating and the 
atoms will scatter off of it. Figure ^| shows the re- 
sulting scattering probability P Pf (t) obtained from this 
approximate model for five atoms and a light field ini- 
tially in a coherent state with mean photon numbers 



N„ 



N-* 



G. 



The atomic Pendellosung oscillations do not decay as 
fast as in the full quantum description of Fig. ITHT a). In 
the present picture, it can be attributed to the degrada- 
tion of the intensity grating. The maximum oscillation 
amplitude occurs when J lies in the equatorial x-y plane, 
but the scattering of the atoms leads to a redistribution 
of the photons between the counter-propagating modes 
and a decrease in the optical fringe visibility. 



5/. 



si 



(St 



2i 



(si 



si 



(15) 



(see Eq. ©). The S l k ,l = x,y,z obey the usual angular 



momentum commutation relations 



\_Sk: S^r^ — itlmn&k,k' S k ■ 



(16) 



where ei mn is the Levi-Civita symbol. Using these com- 
mutation relations we obtain the coupled equations of 
motion for the light field and atomic operators (|10|l : 



IV. STANDING WAVE QUANTIZATION 



For a standing-wave quantization of the light field the 
Hamiltonian Q is replaced by 



dt 



g{S x $L + S v y) x J, 

(6w k & + g(J*Z+jyy))xS k , (17) 



where N — a^a, the number operator for the optical field 
mode, is clearly a constant of motion, and <r and a are 
bosonic creation and annihilation operators. 
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FIG. 10: Scattering probability for five atoms scattering off (a) a running wave and (b) a standing wave in a coherent state. In 
both cases the recoil energy is Ei q = 50g and the Fermi momentum is kF = O.lq. The mean number of photons is N q = TV-, = 6 
in (a) and TV = 12 in (b). Time is in units of <? _1 and momentum in units of q. 



(a) 



(b) 



(c) 



(d) 



0.5 



1.5 2 2.5 

t/27t g~ 1 



3.5 



FIG. 11: Mean number of scattered atoms N ac (t) in coher- 
ent state light fields, (a) Running-wave light field with mean 



photon numbers N q 



N- 



6. The Fermi momentum is 



kF = O.lq and the recoil energy is E2 q = 50g. (b) Light field 
as in (a), atoms without dephasing, i.e. fcj? = 0. (c) Standing- 
wave light field with mean photon number N = 12 and Fermi 
momentum kF = O.lq and recoil energy Ei q = 50<?. (d) Light 
field as in (c), but without dephasing for the atoms, kF = 0. 
Time in units of 2ivg~ 1 . 




FIG. 12: Scattering probability P p , (t) for five atoms in a 
running-wave light field as calculated within the Bloch vector 
picture. The Fermi momentum is kF = 0-lq and the recoil 
energy is E-2 q = 50g. Time in units of <? -1 and momentum in 
units of q. 



A. Raman-Nath regime 



From Eq. I|18f) we now have 
. d 



i—c 
dt 



= (E k2 -E kl )cl iCk2 

gN r ± t \ 

+ — \C kl C k2~2q + C [ kl C k2+2q j 



(19) 



gN 



fei+2g Cfc 2 



Since Eq. i|18fl does not couple states with different 
photon numbers, we can replace N by the correspond- 
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FIG. 13: Scattering of two atoms off of a standing-wave light FIG. 14: Scattering of two atoms off of a standing-wave light 

field in a Fock state with N = 6 photons, the recoil energy field in a coherent state . The mean number f photons is 

is E 2q = and the Fermi momentum is k F = 0.1 9 . Time in giX] the recoil energy is £ 2? = g and the Fermi momentum is 

units of g and momentum m units of photon momentum q. kp = Q lq Time ^ unitg of g -i and momentum in units of 

photon momentum q. 



ing eigenvalue N for a particular number state, \N). We 
can then calculate the evolution of (cj^c?^) for a gen- 
eral state of the field by averaging over the appropriate 
photon number distribution. 

The equations for the first-order moments of the indi- 
vidual atoms are identical to those describing the scat- 
tering of a single atom by a classical light field, if one 
identifies gN/2 with the classical Rabi frequency. This 
follows from the absence of correlations between the light 
field and the atoms, together with the condition q > kp, 
which implies the absence of Pauli blocking. 

The scattering of a single atom by a classical field is 
a well-studied problem. An analytical solution is known 
in the Raman-Nath regime, see e.g. |l8|, l2il |22J. If the 
kinetic energy of the atoms is not negligible, on the other 
hand, one has to rely on a numerical solution. 

Figure ^| shows the results for two atoms in a Fock 
state with N=6 photons. For short times the scattering 
clearly resembles the single-particle behavior, the prob- 
ability of finding an atom in the m-th side mode being 
proportional to J^gNt). Note that this result is iden- 
tical to the approximate first-order calculation for the 
running-wave coherent state of section III. 

The solution for a coherent state is obtained by av- 
eraging over a Poissonian photon number distribution. 
The result is shown in Fig. ^]for a mean photon num- 
ber N — 6, which illustrates the dephasing due to the 
distribution of Rabi frequencies for such a state. 



B. Bragg regime 

As previously discussed, the atoms can now be de- 
scribed in a two-state basis, leading to the effective 



Hamiltonian, see Eq. 110(1 : 

H= £ Ucj k S z k + ^(S l :+S+)\. (20) 

ke[-k F ,kp] K ) 

Evidently, the Hamiltonian decomposes into independent 
single-particle Hamiltonians for the individual atoms, 
H = '}2 k Hk 1 so that the atomic equations of motion 
decouple and can be readily solved analytically. If the 
cavity is in a Fock state, the atoms undergo Pendellosung 
oscillations independently of each other, with a detuning 
given by Slju- We then obtain 

(21) 

as illustrated in Fig. HJb) for five atoms and the light field 
in a Fock state with N = 12 photons. As expected, the 
frequency of the Pendellosung oscillations increases and 
their amplitude decreases with detuning, see Eq. 1(21(1 . 
In contrast to the running-wave case, these oscillations 
remain perfectly sinusoidal for all times. 

The oscillations in the mean number of scattered atoms 
N sc (t) is shown in Fig. E{c) for an atom number of 200. 
All other parameters are as before |2^| . Their decay due 
to the spread in detunings StOk, resembles the dephasing 
of an inhomogeneously broadened ensemble of indepen- 
dent two-level systems, with a dephasing time given by 

27T 

The scattering probability for a coherent state light 
field is shown in Fig. HUf bl and the mean number of scat- 
tered atoms N sc (t) is shown in Fig. Hir e - ). The atomic 
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parameters are the same as for the Fock state calculations 
while the light field has been replaced by a coherent state 
with the mean number of photons N = 12. 

As in the running-wave case, the oscillations collapse 
due to the spread in Rabi frequencies associated with the 
coherent state. In contrast to that former case, though, 
they also undergo a revival after the time T rev i va ] = 
2irg~ 1 , as follows from the fact the difference between the 
Rabi frequencies for N and N + 1 photons is g/2, inde- 
pendently of N. As in the two-photon Jaynes-Cumming 
model, all number states rephase at the same instant 
and the revivals are perfect if the dephasing due to the 
kinetic energy of the atoms can be neglected, as shown 
in Fig lTTT d'). With dephasing included, the revivals are 
only partial, see Fig. Illf cl. 

Another difference between the running-wave and 
standing-wave quantization schemes can be seen in the 
evolution of the scattering probability for the coherent 
state. In the case of running waves the probability of 
finding an atom in the scattered states is about 1/2, in- 
dependently of their detuning, while in the standing wave 
case this probability decreases with increasing detuning. 
This difference underlines the importance of the correla- 
tions between the light field and the atoms. While the 
atoms move independently in a standing wave light field, 
the atoms and the field become an inseparable quantum 
system in the running-wave case. 

V. CONCLUSION AND OUTLOOK 

In this paper we have compared the scattering of ultra- 
cold fermions by quantized light fields composed of "true" 
standing waves and of superpositions of counterpropa- 
gating running waves, both in the Raman-Nath and in 
the Bragg regime. The central difference between the 
two quantization schemes is that the entanglement be- 
tween the light field and the atoms plays a crucial role 
for a running-wave light case, but not for standing-wave 



quantization. 

In the Raman-Nath regime the scattering by a 
standing-wave light field in a Fock state is similar to the 
scattering by a classical field, and can be largely under- 
stood from the results for a single atom. The only many- 
particle effect is a dephasing of the atomic motion due to 
the finite width of the initial momentum distribution of 
the atoms. For running-wave quantization, on the other 
hand, the quantum correlations that develop between the 
light field and atoms are essential, and their neglect leads 
to the absence of atomic diffraction by a Fock state of the 
field. 

In the Bragg regime and for a standing- wave light field, 
atomic diffraction can be solved in terms of solutions 
for single atoms, which undergo Pendellosung oscillations 
that undergo a series of collapse and revivals in the case of 
a coherent light field. For the case of running-wave quan- 
tization, these oscillations decay even for a Fock state, a 
consequence of the correlations that develop between the 
light field and the atoms. A coherent state merely leads 
to a faster collapse. 

Table [I] gives a summary of the mechanisms that lead 
to collapses and possibly revivals in the various cases that 
we have disccussed, as well as the associated time scales. 

This paper has only considered fermionic atoms. 
Bosonic operators commute rather than anticommute, re- 
sulting in different equations of motion for the second- 
and higher correlation functions of atomic operators. 
Consequently, differences in the scattering result solely 
from the effect of these higher-order correlations on the 
dynamics. 
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